Two-channel Kondo physics in two-impurity Kondo models 
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We consider the non-Fermi liquid quantum critical state of the spin-S two-impurity Kondo model, 
and its potential realization in a quantum dot device. Using conformal field theory (CFT) and the 
numerical renormalization group (NRG), we show the critical point to be identical to that of the 
two-channel Kondo model with additional potential scattering, for any spin-S*. Distinct conductance 
signatures are shown to arise as a function of device asymmetry; with the 'smoking gun' square-root 
behavior, commonly believed to arise at low-energies, dominant only in certain regimes. 
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Systems comprising several quantum impurities inher- 
ently display an interplay between impurity-bath and in- 
terimpurity couplings [1-7]: the tendency for impurities 
to be Kondo-screened by conduction electrons competes 
with screening by interimpurity spin-singlet formation. 
Rich physics can thereby arise, as demonstrated for ex- 
ample in coupled quantum dots 0, [1| > magnetic impuri- 
ties in metals and recent two- impurity STM exper- 
iments Indeed, the same essential physics governs 
the analogous propensities for heavy Fermion behavior 
or magnetic ordering in lattice systems (lo| . 

The two-impurity, spin-^ Kondo model (2IKM) is the 
simplest to capture this competition local singlet for- 
mation is favored by an interimpurity exchange while 
coupling of each impurity to its own conduction channel 
favors separate Kondo screening below an effective single- 
channel, single-impurity scale Tk 0-13 ■ The lack of inter- 
channel charge transfer in the 2IKM permits two distinct 
phases, and a quantum phase transition (QPT) results 
on tuning Kc ~ Tk- At the critical point, non- Fermi liq- 
uid (NFL) physics arises below Tc, characterized [3] by 
anomalous properties such as fractional residual entropy 
and singular magnetic susceptibility. 

This critical physics is also surprisingly robust to some 
perturbations, notably breaking of mirror (parity) sym- 
metry or particle- hole symmetry Such perturbations 
are marginally irrelevant in the sense that the interimpu- 
rity coupling can be retuned to recover the critical point 
in all cases. But despite considerable effort, 2IKM critical 
physics has proved experimentally elusive — mainly due 
to interchannel charge transfer which smooths the QPT 
into a crossover Regular Fermi liquid (FL) physics 
then sets in below an energy scale T* ; and if the degree of 
charge transfer is large enough that T* ^ Tc, no evidence 



of the critical point will be observed . This is the sit- 
uation relevant to the recent two-impurity experiments 
of Ref. 0: coupling between one impurity on a metal 
surface and one on an STM tip was also accompanied by 
strong tip-surface tunneling. 

Reducing the degree of interchannel charge transfer 
might be possible in a quantum dot device such as that 



proposed in Ref. [i]. Provided T* < T^, NFL behavior 
should be observable in an intermediate energy window, 
as can be understood from a 2IKM critical perspective 
(indeed the eventual crossover to FL physics is wholly 
characteristic of the intermediate NFL state (l2|). 

An alternative route could however involve use of a 
quantum box, which acts as an interacting lead 
Coupling a single dot to one regular lead and one box 
tuned to the Coulomb blockade regime, suppresses in- 
terchannel charge transfer completely. This has been 
exploited to access single-impurity two- channel Kondo 
(2CK) physics [H, 15| in a real device Here we pro- 
pose simply to interject a second dot in series between 
the 'leads' to realize 2IKM physics. While parity and 
particle-hole symmetries are thereby broken, the QPT 
itself is unaffected. Robust NFL behavior should persist 
down to the lowest energy scales at the critical point. 

Here we address two key questions in regard to poten- 
tial realization of 2IKM physics. First, what is the nature 
of the critical point itself? We show that it is identical 
to that arising in a 2CK model with additional potential 
scattering, independent of parity breaking. Further, we 
show that the same QPT and 2CK critical point arises 
in the spin-S" generalization of the 2IKM. Second, what 
are the signatures of criticality in measurable quantities 
such as conductance? These reflect RG flow from higher- 
energy fixed points (FPs) , and depend sensitively on par- 
ity breaking. We find in particular that the square-root 
behavior commonly anticipated 0] at low-energies, is 
absent in the standard channel-symmetric 2IKM. 

Nature of 2IKM critical point.- The 2IKM reads: 

H21K = Ho+Hps+JlSl-sol + JrSb-sqr+KSl-Sr, (1) 

where Hq = J2a.k '^k'^kaa'^kaa dcscribcs two free conduc- 
tion channels a — L/R, with density of states p, and spin 
density at the impurities soa = J2aa' ^laa{\^'y'^')'^oa'a 
(where V^J^tq = ^k^Aaa)- Potential scattering is in- 
cluded via Hps = La ^"V'oaaV'offQi ^ud Sa are spin-i 
operators for the impurities. The 2IKM has been exten- 
sively studied using a number of powerful techniques 0- 
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Figure 1: Spectra Dpa{ijj) vs frequency |i^|/-D for channels 
a — L and R [upper and lower panels] for the asymmetric 
2IKM with fixed Kondo exchanges pj^ = 0.1, pJr ~ 0.05, 
varying interimpurity exchange K (and Hps = 0). Plotted 
using K — Kc{l ± 10~") for dashed and solid lines, with 
integer n = 1 — > 6 approaching progressively the critical point 
^Tk ^ 10-^ D (dotted line). Circles: for a ICK model 
with pj = pJl- Diamonds: Dp2CK{^) for a pure 2CK model 
with Tjf^ = Tc. Squares: Dp2ck{uj) = 1 - Dp2ck{lo). 




and certain similarities have been found 0- 
id l between it and the 2CK model 14, 15| . 



H2CK — Hf) + Hps + JlS ■ Sol + JbS ■ sqr, (2) 



with S a spin-i operator for a single impurity, exchange- 
coupled to two independent conduction channels. The 
physics of the 2CK model is itself immensely rich [l^ : the 
impurity is fully Kondo screened by the more strongly- 
coupled conduction channel below an effective single- 
channel scale Tk, producing two distinct phases as a 
function of {Jl — Jr)- But when Jl = Jr, the frus- 
tration inherent when two channels compete to screen 
the impurity results in 'overscreening', and NFL physics 
results below T^P^ Strikingly, both 2CK and 2IKM 
have the same fractional residual entropy, S'imp = ^ lii(2). 

CFT has been used to describe the critical points of 
both models 0, ■ In the 'unfolded' representation, Hq 
is written in terms of left-moving chiral Dirac Fermions. 
The CFTs for each model can be separated into dif- 
ferent symmetry sectors. In particular, the SU{2)2 x 
SU {2)2 X U{\) flavor, spin and charge symmetries of the 
2CK model [ill, and the ;7(1) x U{1) x SU{2)2 x Z2 
left /right charge, total spin and Ising symmetries of the 
2IKM [il can be exploited. The 2CK and 2IKM criti- 
cal fixed point Hamiltonians take the same form as iJg, 
but with modified boundary conditions (BCs) that affect 
only the spin sector of the 2CK model [H| or the Ising 
sector of the 2IKM The finite size spectrum (FSS) 
at the critical point of each model, can then be deter- 



ps — 0, the FSS of the 2CK critical point is 

I, 1...; while for 



and for 

characterized by the fractions 0, i, i, 
2IKM a different FSS arises: 0, |, i, |,"l... 

Despite these apparent differences between the criti- 
cal FPs of the two models, the 2IKM can be mapped 
onto an effective 2CK model in special cases [1, H, [ll| . 
The key requirement for that mapping is of course the 
generation of an effective spin-i local moment (LM), 
which can then be overscreened by symmetric coupling 
to two conduction channels. In the channel-asymmetric 
limit Jl ^ Jr, 2CK critical physics arises via a simple 
mechanism first involving Kondo screening of the L- 
impurity by the L-lead on the single-channel scale T^; 
followed by second-stage overscreening of the i?-impurity 
by the i?-lead and an effective coupling to the remain- 
ing Fermi liquid bath states of the L-lead below Tc- An 
effective 2CK model of form Eq. [51 valid at low-energies 
T < Th, can be derived formally using the approach of 
Ref. [l3| , exploiting the Wilson chain representation [l8| ; 
and effective couplings follow as pJ^ ^ K/T^ and 
pJr ^ [^/pJr - ^IpJlY^- The 2CK FP is thus sta- 
ble when K = Kc ^ T^pJ"^, so the low-energy physics 
of the 2IKM is in this case wholly equivalent to that of 
the 2CK model. 

Importantly however, the L-channel free electrons in 
the effective 2CK model acquire a 7r/2 phase shift due 
to the first-stage single- channel Kondo screening of the 
L-impurity in the original 2IKM. This is seen clearly in 
the dynamics of the asymmetric 2IKM; to demonstrate 
which, and to highlight the basic physical picture. Fig. [1] 
shows spectra Dpa{uj) = — 7rpIm[tQ(w)] vs |w|/Z3 with 
ta{ijj) the scattering t matrix \\^. Results are obtained 
from NRG, exploiting all model symmetries, discretizing 
conduction bands of width 2D logarithmically using A = 
3, and retaining 8000 states per iteration in each of z = 3 
interleaved calculations (for a review, see Ref. [isj). 

Single- channel Kondo screening of the L-impurity by 
the L lead on the scale of is seen directly in the L 
spectra in the upper panel: a Kondo resonance, reach- 
ing the unitarity limit DpL = 1, which has precisely the 
form of a regular single-channel Kondo (ICK) model 
This embodies the 7r/2 phase shift in the L channel; but 
no such feature is observed on this energy scale in the 
R spectra (lower panel), indicating that the i?- impurity 
is still essentially free. On tuning the interimpurity cou- 
pling K closer to the critical point of the 2IKM, the spec- 
tra in both channels fold progressively onto the critical 
spectra. For energies \w\/D <C T^, the critical spectrum 
Dpr{uj) is precisely that of a 2CK model Dp2CK{^) with 
— Tc- But to leading order the critical s pec trum in 
the left channel is Dpl{uj) = 1 — Dp2CK^) [III- Thus, 
at the channel-asymmetric critical point 



Dpa{uj) ~ 2 



(3) 
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with a = ±1 for channel L/R resulting from the addi- 
tional L-channel 7r/2 phase shift, and (5 a constant 0(1). 

This phase shift can be included in the 2CK model, 
Eq. [51 via the potential scattering term Hpg (ie Vl ^ oo 
but Vr = 0, accompanied also by retuning Jl and Jji to 
access the critical point). This is equivalent to adding in- 
finite uniform and staggered potential scatterings, which 
affect respectively the charge and flavor sectors of the 
2CK model. Modifying the CFT for the critical point 
of the 2CK model to include this, we find [20] the BC 
becomes equivalent to that of the 2IKM. The FSS is also 
naturally affected and is given by (20| 



E. 



2CK 



\{Q - af + jJU + 1) + \jf{]f + 1) - hfp, (4) 



where Q, j and jp are the charge, spin and flavor quan- 
tum numbers. Uniform potential scattering shifts the 
charge parabolas, while staggered potential scattering 
biases the flavor sector. The 7r/2 phase shift in the L- 
channel corresponds to a = 1 and b — ^ [20]. Only cer- 
tain quantum number combinations are allowed at the 
critical point, as given by the nontrivial gluing condi- 
tions derived in Ref. [31; and which reproduce fully the 
2IKM spectrum when used with Eq.|4](see Fig. 2 of [l^). 
One remarkable result obtained from our NRG calcula- 
tions ji^] is that the FSS at the critical point of the 
2IKM does not depend on channel asymmetry (whence in 
particular the critical point possesses an emergent parity 
symmetry, irrespective of bare model symmetries). Fur- 
ther, we have shown [20] that the critical point for one 
model with potential scattering Vl and Vr is equivalent 
to the critical point of the other model with different po- 
tential scattering Vl and Vr . The 2IKM and 2CK critical 
FPs are thus equivalent in the sense that they lie on the 
same marginal NFL manifold parametrized by Hps . 

Conductance lineshapes and symmetry. - Full RG flow 
from the local moment (LM) FP to the 2CK FP is 
thus recovered at the critical point of the asymmet- 
ric 2IKM. This is manifest [6] in the conductance aris- 
ing e.g. when a given channel a = L/R is split into 
source and drain. At zero-bias, it is given exactly [21 ] 
in terms of the scattering t matrix (considered for the 
channel-asymmetric case in Fig. [1]) by G2jj^{Vsd = 
0,T)/(2e2/^-iGJJ) = - J^^dudfioj/T)/dL,Dp^{u,T); 
with f(uj/T) the Fermi function, and the impurity-lead 
coupling parametrized by = 4r^rjJ/(r^ + T^)^, in 
terms of the a = L/R hybridizations to source (F") 
and drain (F^). Indeed, in the limit F^ > (ie. 
Gq ^ 1), the T = conductance follows as G2//f (V^sd) = 
G2iK(ysd,T = 0)/i2e^h-^G^) = Dp^{u = Vsd,T = 0), 
and hence from Eq. [3]one finds at low energies Vsd ^ 7c, 



G'^iKiVsd) = \ + a/3v/VW7; + 7„(^sd/Tc) + ... (5) 

where we include also a term linear in Vsd/Tc- 

The leading square-root behavior of Eq. [S] has been 

viewed as the 'smoking gun' signature of this 2CK 
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Figure 2: T = conductance G21K / i'^e^h-^Go) through 
channel ct = L and R (solid and dashed lines) vs bias Vsd/ D, 
at the critical point. Shown for pJl = 0.075 > pJn, varying 
pJr = 0.075 0.05 in steps of 0.0025, with K = Kc T^i^ 
retuned in each case (and H^s = 0). Thick solid line is the 
symmetric case Jl = Jr\ asymmetry Jl/Jr > 1 increases 
in direction of arrows. Circle and diamonds: pure ICK and 
2CK scaling spectra. Inset: zero-bias conductance vs T/Tc 
for Jl/Jr = 1 and 2 (and a = L, R), exhibiting respectively 
leading linear and square-root behavior (dotted hues). 



physics (1, [13, and was used to identify the critical point 
in the 2GK experiment of Ref. [13] . But we note that, 
unlike the 2CK model, the 2IKM does not possess SU{2) 
flavor symmetry. Since symmetry dictates which opera- 
tors can act in the vicinity of the critical FP, this is natu- 
rally reflected in the asymptotic conductance through the 
coefficients /3 and 7^. Indeed, the full energy-dependence 
of conductance depends on the unstable FPs, whose vy- 
ing effects on RG flow again depend on symmetry and 
model parameters. For example, in the usual symmetric 
2IKM (Eq. [T] with Jl — Jr and Hps = 0) , no incipient 
LM is formed: there is no intermediate energy window 
with eg. Simp = lii(2) entropy, and RG flow proceeds 
directly to the 2CK FP from the LMxLM high energy 
FP describing a pair of free impurities (Simp = ln(4)). 

The effect of parity-breaking is explored in Fig. [2l 
showing NRG results for conductance vs bias Vsd at 
the 2IKM critical point (obtained for Gq <^ 1, as 
above). Conductance in the asymmetric limit Jl S> Jr 
is consistent with Ref. [^, and physical expectation as 
above (see Eq. [5]). Here, the asymptotic conductance is 
G2iK(^sd) ^ l-G^jK{Vsd) = G2CK{Vsd) at low energies 
Vsd « T^, where T^:$>Tc = T]f^ and with Gzc^V.d) 
the conductance of the standard 2CK model [l3l.l2^ (di- 
amonds). Thus, on exchanging Jr ^ Jl, the coeffi- 
cient f3 of Eq. [5] must change sign. But what happens 
as the asymmetry is decreased? We find the leading 
square-root contribution in Eq. [S] vanishes (see Fig. [5]) , 
as /3 {Jl ~ Jr), and leading linear behavior emerges at 
the symmetric point Jl = Jr (the same naturally arising 
as a function of T at zero-bias, see Fig. [5] inset). In fact 
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Figure 3: Impurity entropy S'imp(r) vs T/D for tlie 2IKM 
with spin-1 impurities. Plotted for fixed pJh = 0.15, pJn = 
0.05, varying K/D = 0, 10"^^ 3 x 10"'', ~ 6 x 10"^ 2 x 
10~* and 10~^ for lines (a)-(f). The corresponding physical 
processes are illustrated in panels (a)-(f): impurities denoted 
as circles and conduction band orbitals (in the Wilson chain 
representation [l^ ) as squares. For discussion, see text. 



Vmeai-Vsd behavior also emerges as the symmetric point 
is approached, since the square-root term dominates over 
a shrinking window Vsd/Tc <C (/3/7q)^. 

Another striking feature of the conductance in more 
channel-symmetric situations is the behavior at higher T 
or energies > Tc w T^. Here the behavior is wholly 
characteristic of single-impurity, single-channel Kondo 
physics, as seen by comparison to the circles in Fig. [21 

The absence of square-root behavior in conductance of 
the symmetric 2IKM is contrary to common belief 0| , 
so we sketch now our CFT proof As pointed out in 
Refs. 0, corrections to the t matrix in the vicinity of 
the critical point (whose oj-dependence displays the same 
scaling as conductance) are determined from irrelevant 
boundary operators consistent with symmetry. Two such 
play a role here: 5Hi = cie' and SH2 = C2 J-i • (in the 
notation of Ref. Q). The operator jli • is the leading 
irrelevant operator of the 2CK FP, whose effect on the 
t matrix is known [2^ to yield the famous square-root 
behavior. However, J-i ■ (p has odd parity in the 2IKM 
(unlike 2CK), which implies that its coefficient C2 ~ ( Jl — 
Jr) vanishes in the symmetric limit. In both models, 
SHi does still contribute (ci always being finite 01 ). One 
might naively expect e' to behave similarly to J-i-4> since 
they have the same scaling dimension 3/2. However, the 
key difference between e' and J_i -0 is that only the 
latter is a Virasoro primary field. In consequence [20|, the 
leading square-root correction to the t matrix from 5Hi 
vanishes. In the symmetric 2IKM, the leading square- 
root behavior of conductance thus also vanishes. 

Spin-S 2IKM.- Multilevel quantum dots can behave 
like S = 1 impurities (23 |. and high-spin impurities such 
as Co {S = 3/2) have been manipulated with STM 



Thus a natural and pertinent generalization of the 2IKM 
involves spin-S* impurities: the model remains Eq. [1] but 
Sr,Sl are now spin- 5 operators. 

A QPT must again arise, as follows from the same 
line of argument as the spin-i 2IKM On tuning K 
there is a phase-shift discontinuity on going from the lo- 
cal singlet phase for large K to a separated spin-S" un- 
derscreened Kondo phase for small K . The nature of the 
transition arising at Kc is again clear by considering the 
asymmetric limit Jl 3> Jr. For concreteness consider 
S = 1, although the argument extends easily to higher- 
S. NRG results for the entropy S'imp(r) vs T are shown 
in Fig. [3l together with cartoons highlighting the key 
physical processes. In (a) the impurities are completely 
decoupled {K = 0), with each thus underscreened to a 
spin-i by its attached lead a, on its own single-channel 
Kondo scale (with residual entropy 2 ln(2)). For small 
finite K < Tj^, (b), these residual moments form a local 
singlet state on the scale T ^ K, so the residual entropy 
is quenched. On increasing the interimpurity K further 
{T§ < K < T^), the underscreened spin-i L impurity 
and the still unscreened spin-1 R impurity are coupled 
and form a local doublet state on the scale T ^ K. This 
can then be single-channel Kondo screened by an effective 
coupling either to the L channel (c) or the R channel (e) , 
and the residual entropy is again quenched. However, 
L and R effective couplings can become equal on fine- 
tuning K. This is the single spin-i 2CK critical point 
(d), with residual entropy iln(2). For large K ^ T^, a 
local interimpurity singlet state arises as expected (f). 

Analysis of the finite size spectrum at the critical point 
shows it to be identical to that of the regular spin-^ 
2IKM, independent of asymmetry and is hence that 
of a 2CK model with additional potential scattering. 

Conclusion.- We have shown the critical point of the 
spin-S 2IKM, including the spin-i variant, to be ubiqui- 
tously 2CK in nature. However conductance lineshapes 
measurable in experiment exhibit distinctive behavior de- 
pending on underlying symmetries, the low-energy be- 
havior in particular evolving from square-root to linear 
behavior in Vsd or T as the channel symmetric point is 
approached, and for any spin-S*. 
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I ! 
O 

^ ■ I. CONNECTING THE FINITE SIZE SPECTRA OF THE TWO-CHANNEL AND TWO-IMPURITY 
^ ; KONDO MODELS 

o . 

, In this appendix we demonstrate that the critical points of the 2CK model and 2IKM are connected by a marginal 
operator corresponding to the potential scattering Hpg . Specifically, the finite size spectrum of one model is shown to 
[ transform into that of the other model when this additional term is accounted for explicitly. In the channel-asymmetric 
^ ' 2IKM, the physical origin of this operator is readily understood from the mapping^ between 2IKM and 2CK, where 
an additional phase shift in the more strongly-coupled channel arises due to first-stage single-channel Kondo screening 
. of one of the impurities. 

First, we review how the finite size spectrum at the critical point of the 2CK model is organized according to the 
unperturbed CFT with [/(I) x SU {2)2 x SU {2)2 charge, spin and flavor symmetry sectors. The spectrum consists 
^ • of Kac-Moody conformal towers, obtained by combining these symmetry sectors according to the gluing conditions^ 



O 
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shown in Table ID The energy of the lowest lying state in each conformal tower (measured in units of ^^^"-^ , with 
vp = \ the Fermi velocity and ft = 1 hereafter) follows from the formula^ 



pH . i^ICK — — H ] 1 ^ , 

o 

. where Q is the total charge, j = 0, i, 1 is the total spin, and jp = 0, i, 1 is the flavor quantum number. Only certain 
\ quantum number combinations are allowed at the 2CK critical point^ as given by the nontrivial gluing conditions of 
I ■ Table H] (and with level degeneracies also following from the Table). Also note that the gluing conditions for integer Q 
K*" ' are defined modulo 2. Each conformal tower has an infinite number of states. Consider for example the spin j tower of 
. the spin SU {2)2 sector. The states with lowest energy form a spin j representation of SU(2), and hence are {2j + 1)- 
degenerate (the corresponding magnetic quantum number is jz = —j, —j + 1, The reducible representations of 

SU{2) with increasing spin (equal to j modulo 1) give rise to higher energy states. In the space of jz and energy n, 
the envelope of an SU{2)k conformal tower on level k with spin j is n = [{j^Y ~ J^]/^- For ^ visualization of such 
unperturbed conformal towers, see the upper panels of Fig. [TJ 

We now perturb the CFT by introducing potential scattering Hps = ^a'tpoaa'^oaa- To see how this perturbation 
enters in the CFT, we first switch off the Kondo interactions, and consider free Fermions of a single species aa. With 
antiperiodic boundary conditions ipcyai—L/2) = —ip^aiL/2), the single particle momenta are k = ^{n + i); and 
for linear dispersion, the single particle energies are just given by k. The potential scattering amplitudes Vq. can be 
parametrized in terms of scattering phase shifts 6a- They are defined in terms of the shift of these single particle 
$—1 ' levels k = + ^ — For A^^a particles of each species, the correction to the total energy (again in units of ^) 
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TABLE I: Nontrivial gluing conditions of conformal towers belonging to the U(l) x 5(7(2)2 x 517(2)2 charge, spin and flavor 
symmetry sectors. 



FIG. 1: Top panels: SU{2)k Kac-Moody conformal towers for k — 2, for j — 0,^,1 from left to right. Bottom panels: tilted 
towers due to Zeeman term —bjp with 6 = |. Multiplicities are not shown. 



is thus 
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In terms of the particle numbers, the total charge is Q = J^ai-^o'L + ^o-fl) while the magnetic quantum number of the 
flavor sector is jp = \ '^^{N^l — N^r). In the presence of Hps, the spectrum thus follows as E2CK ~ E2CK + AE, 
viz 

E,cK - ^ + + ^ - 

where S^a = Sa since the potential scattering Hps is independent of spin a. Eq. ^ is equivalent to Eq. (4) of Ref. [l|, 
with a — 2^^^i^ and b = ■ Thus Hps simply biases the charge and flavor conformal towers. For a visualization 

of the effect on the SU{2)2 flavor towers, see the lower panels of Fig. [1] The spectrum obtained using the trivial 
free Fcrmion gluing conditions^ describes free electrons but with the additional potential scattering. Importantly, 
switching on Hps and switching on the Kondo interactions commute, since there is a spin-charge-flavor separation in 
the modeli^ Employing instead the nontrivial gluing conditions of Tabic |T] in Eq. ([3]), we thus obtain the 2CK critical 
spectrum including the influence of Hps- 

In Fig. [21 we examine the evolution of the energy levels at the critical point of the 2CK model as function of 6l 
for Sr = 0, plotting the lowest excitation energies 6Ei ^ Ei — min{Ei}. The fractions 0, i, i, |, 1... correspond to 
multiplets of the symmetry U{1) x SU{2)2 x SU{2)2 of the regular 2CK model. Those multiplets split due to the 
flavor field and charge field induced by finite Hps, but eventually at 6l = 7r/2 they recombine to give new multiplets 
characterized by the fractions 0, |, i, |, 1..., and which correspond to that of the regular 2IKM<^ 

To strengthen the connection between the two models, we now consider the effect of Hps on the critical CFT 
of the 2IKM itself. We exploit here the Bose-Ising representation^ of the 2IKM, consisting of a decomposition into 
SU{2)i X SU{2)i X SU (2)2 X Z2 left/right charge (isospin), total spin and Z2 Ising symmetry sectors. The ^-component 
of a = L or i? isospin is essentially the total charge of channel a. Specifically, its magnetic quantum number is given 



5 ^ in ^ 

8 4 T 2 

Sl 

FIG. 2: Excitation energies 5E according to Eq. ((3| as a function of 5l for 5r = 0, connecting the regular 2CK spectrum {5l = 0) 
and the 2IKM spectrum {5l ~ 7r/2). Note that at Sl ~ tt/8 the ground state changes from the doublet {Q,j,jF) = (0, 1/2,0) 
to the jp = 1/2, (5 = 1 component of the quintet (±1, 0, 1/2). 



simply hj ia = \ X^cr ^<ya. in terms of the particle numbers. Using Eq. ©J it follows that /S.E = — 5aia- The 
energy of the lowest states of each conformal tower in the 2IKM with the additional potential scattering term Hps is 
thus E21K = E21K + where E21K is given in Ref. Our final result is 

~ _ »l(»l + 2) iRjiR + 2) + Jl , Jr 

E2IK - 5 1 5 \ h Xlsing - 2 - 2 1^, (4j 

where, = 0, ^ are the isospin quantum numbers, j = 0, ^, 1 is the total spin quantum number and xising = 0, j^, ^ 
(corresponding to the unity operator, spin field a, and fermion field e of the Z2 Ising model). The nontrivial gluing 
condition of the SU{2)i x SU{2)i x 51/(2)2 x Z2 sectors is given in Table. Ill of Ref. and gives the finite size 
spectrum of the 2IKM in the presence of potential scattering. 

Using Eqs. ([3]) and (|4]) with their corresponding gluing conditions, we find 

n'S'iSL,6R) = nl^^'iSL + 7r/2, Sr), (5) 

where 'Hcp = {^Ei} at the critical point of either model, parametrized in terms of the free Fermion phase shifts Sl 
and Sr. The finite size spectrum at the critical point of one model with potential scatterings Vl and Vr is identical 
to that of the other model with different Vl and Vr when the phase shifts satisfy Eq. ([5]). Thus, the critical FPs are 
equivalent in the sense that they lie on the same marginal manifold parametrized by Hps. 

A. Alternative formulation in terms of Majorana Fermions 

Having drawn the connection between the non-abelian CFT formulation of the 2CK and 2IKM, we discuss now the 
same relation within the framework of abelian bosonization. 

The steps applied first by Emery and Kivelsor^ as an alternative solution to the 2CK model, and later by Gan^ 
in the 2IKM context, are as follows, (i) Bosonize separately each chiral Fermionic species (with the impurities at 
the 'boundary' located a.t x = 0), ^/jaa{x) ^ g-i4>ac(x) ^ ^jj^ Define linear combinations of the bosonic fields (j)aa{x), 
{0c, 0s, 0/, 0x} = \ Z)^,Q=i 4><ya{l, (-l)'^+\ (-l)"+\ (-1)''+"}, corresponding to charge, spin, flavor and spin-flavor 
bosons {A = c,s,f,X). (iii) Define four new Fermionic species by ^Aix) ^ e"*"^-*^^-*, whose real {xt = ^^14+^) and 
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imaginary {xf = ) parts give eight Majorana Fermions (MFs) xf{^)j (with j = 1, 2 and A = c, s, /, X), which 

may be regarded as components of the vector xi^)- The Hamihonian for the free theory with = Ji?, = and 
Hps = is then H = ^ dx x{^) ' dxx{x)^ with the trivial boundary condition (BC) x(0~'') = x(0~)- The BC 
relates the field before scattering at a: = 0+ and after scattering at a; = 0^ (using a left moving convention). The 
remarkable fact is that in both 2CK and 2IKM, the critical FP Hamiltonian is the same as that of the free theory, 
but with a modified BC that is also simple in terms of the MFs. For 2CK it is simply^ 

Xs(0+) = -x.s(0-) and x/(0+) - x/(0~) for (A, j) ^ (s, 1), (s, 2), (X, 1), (6) 

with Xs = {XiiX2iXi)i such that the modified BC affects only the spin sector of the 2CK model. For 2IKM)^ 

X^(0+) = -X^(O-) and x/(0+) = xf (0") for (Aj)^(^,2), (7) 

where the modified BC here affects only the Ising sector of the 2IKM. 

An odd number of MFs suffer the modified BC in both models. This is a hallmark of NFL behavior, as it implies^ 
a vanishing amplitude for scattering of an electron into an electron. By contrast, the BCs for Fermi liquids always 
involve an even number of MFs, corresponding thereby to linear BCs for regular Fermions. 

Besides their common NFL character, the different BCs for the critical points of the 2CK and 2IKM do not 
immediately suggest any connection between the models. However, we now show that adding potential scattering in 
one channel of the 2CK model (Eq. (2) of Ref. [ll, with Vl ^ oo and Vr — > in the added term Hps) causes a change 
in the critical 2CK BC which makes it equivalent to that of the 2IKM. 

As above, the addition of infinite potential scattering in the left channel (which results in the desired 5^ = 7r/2 
phase shift) is equivalent to the addition of infinite uniform potential scattering + ■^ipR'^'R = ^xiXi which 

affects only the charge sector, and infinite staggered potential scattering ^ipl^ipi^ — ^V'ijV'i? = ixixij affecting only 
the flavor sector. Since both charge and flavor sectors are unaffected by the Kondo interaction (which acts purely in 
the spin sector) , the resulting behavior can be understood from the response of free Fermions to potential scattering 
in the simpler situation Jl = Jr = 0. For — >■ oo and Vr 0, modifled BCs are conferred to the charge and flavor 
MFs: x?{(0+) = -Xii{0~)- Together with the BC from Eq.H we now obtain the BC for the critical 2CK model in 
the presence of the potential scattering, viz 

Xf(0+) = xf(0-) and x/(0+) = -x/(0-) for {A,j)^iX,2). (8) 

This BC is the same as that of the 2IKM, Eq. (O, except for an additional minus sign afflicting the BC for all MFs. 
However, the BCs corresponding to Eqs. ([7]) and ([H]) are in fact equivalent in the 2CK model, as now shown. 

An 'overall' minus sign in the BC for the MFs can be interpreted as a phase shift felt by each of the charge, spin, 
flavor and spin-flavor Fermions, deflned in terms of the scattering process V'a(O^) = cxp(2iSA)'4'A(0~^), and with 
Sa ~ 7i'/2 for A ~ c,s,f,X. Since the number operators for ipA Fermions are related^ to the number operators 
for the regular V^a Fermions via {Nc,Ns,Nf,Nx} = i E<Ta ^'^"{l' (-l)''"^^ (-l)""^^ (-1)""^"}' follows that 
'^j^Na = 2N^L. In terms of the particle numbers Na, the correction to the total energy due to additional phase 
shifts is AE = — i J2a ^aNa — ~N^l- Comparing to Eq. [21 we then find that J-fL = tt (and Jo-q = for era L). 
Thus the additional global minus sign in the BCs for the MFs equates to V'ti(0 ) = exp(2i(5-^i)?/'tL(0^) = ij^^hi'^^) 
for the original Fermions, and so has no effect on the finite size spectrum — or on the lowest-energy physics in general. 

II. CORRECTION TO THE GREEN FUNCTION NEAR THE 2IKM CRITICAL POINT 

In this appendix we consider corrections to the critical FP Hamiltonian of the 2IKM, and from them construct 
corrections to the Green function. We show that the nature of these corrections (which show up in physical quantities 
such as conductance) is characteristic of the NFL behavior, and depend distinctively on parity breaking in the model. 

A. Fixed point Hamiltonian of the 2IKM 

First we construct the list of operators allowed by symmetry in the generic asymmetric 2IKM (see Eq. (1) of Ref.[l|), 
including the dependence on small Jl — Jr and on small potential scattering amplitudes Vl and Vr. We rely on the 
same operator content obtained by Affleck and Ludwig;^ based on the SU{2)i x SU{2)i x SU{2)2 x Z2 symmetry 
decomposition. The critical FP of the symmetric 2IKM without potential scattering has only one irrelevant operator 
of dimension 3/2, denoted e', and is given by Eq. (5.1) of Ref. [H. In the present context where parity and particle-hole 
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TABLE II: Operator content of the 2IKM critical FP. 

symmetries are in general broken, additional operators are allowed, and in fact dominate the approach to the NFL 
FP. 

Each operator that can occur at an unstable FP is either a primary operator or a descendant of it.— The primary 
operators are labeled with the same quantum numbers (zi, 12, j, Ising) used to label states in Eq. The particular 
combinations of SU(2)\ x S'C/(2)i x SU{2)2 x Z2 left/right isospin, total spin and Ising operators, together with 
their scaling dimension x, was obtained using the double fusion method in Ref. HI. For clarity and completeness we 
reproduce that list here in Table HTl 

In general, the presence of a given symmetry implies that the only allowed operators are those which are singlets 
of that symmetry. In particular, the full SU{2)i x SU{2)i x SU{2)2 x Z2 symmetry of the standard 2IKM with 
Hps = gives rise to the leading singlet operator (ii, Z2, j, Ising) = (0, 0,0, e). This relevant operator is associated^ 
with the perturbation {K — Kc), which destabilizes the 2IKM critical point. At the critical point K — K^, one leading 
irrelevant operator— is the descendant of (0,0, 0,e), denoted Ci = e' « dre- Since e has scaling dimension 1/2, its 
first descendant has scaling dimension 3/2. But singlets can also be obtained as descendants of primary SU{2) fields 
with integer spin. Indeed, the other singlet operator with dimension 3/2 is in this case obtained by acting on the 
4> = (0, 0, 1, 1) primary vector field of the j = 1 spin sector with an operator that creates spin excitations.— Thus, 
one must consider also a second irrelevant operator O2 = J-i ■ (where J_i is the lowest Fourier mode of the spin 
current creating such excitations). 

Since J_i • (f) is also the leading irrelevant operator in the 2CK model,— one might expect the physical behavior 
stemming from it to be common to both 2CK and 2IKM. However, the crucial difference between the models is that 
J_i • has odd parity in the 2IKM but even parity in the 2CK model; as now shown. 

Parity symmetry (corresponding to left/right reflection in space) is defined by invariance under the permutations 
Sl ^ Sr and ■i/'L i^R- One rather subtle consequence^ of parity symmetry in the 2IKM is that 4> —4>. This 
can be seen from the parity-odd operator ^\dipj^ — iPr<^'4'b, ~ (and noting that e is the operator corresponding 
to {K — Kc), and is as such unaffected by the reflection, e o e). Thus J_i • (f> is also parity-odd in the 2IKM. But 

in the 2CK model, ipl.'^'^L ~ '^^r^'^r ~ $4''f- Hsre, $f is the jp = 1 primary field of the flavor sector, and <p^p is its 
z-component. Since parity corresponds to a tt rotation around the x flavor direction in the 2CK model, the flavor 
current Jp = ^pj^^ip^ transforms as Jf. — > Jf., Jp^ —Jp'^', and similarly 0^ — > (j)p, cjy^p^ — > —(j)^p^. Thus J_i • (j) 

is parity even in the 2CK model. Of course, this is the expected result here since J_i and are both spin operators, 
but the parity transformation affects only the flavor sector of the 2CK model. 

Finally, we consider the additional operators appearing when isospin SU{2)i x SU(2)i symmetry is broken (but 
spin SU(2) symmetry is maintained), as occurs when potential scattering from Hps is included. Specifically, Hps = 
2 Val^{x = 0), where the z-component of the channel-a isospin current is given by I^i^) = ^ X^cr '^Xa{^)'4'(ya (^)- 
Since Va does not couple to fields with half-integer isospin, and the primary operators have total spin j = 0, then it 
follows that Va couples only to dimension 3/2 descendants of the (0, 0,0, e) field. Thus, two further dimension 3/2 
operators ©3,4 — {II fl)-ie are obtained (where (/q)-i is the isospin analogue of the spin operator J_i discussed 
above). No other dimension 3/2 operators arc consistent with the symmetry. 
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The corrections to the critical FP Hamiltonian in the presence of Hps thus comprise four operators, 

As shown in Ref. the coefficient Ci ~ 1. By contrast, the odd transformation property of O2 under parity imphes 

for small small ( Jl — Jr)- Finally oc Vl,r for small Vl,r- 

In the conformal limit (which relies upon linear dispersion and infinite conduction bandwidth), the Bosc-Ising 

representation of the 2IKM reads H21K = Hq + J(0) • {JlSl + JrSr) + 0(O)e • (JlSl ~ JrSr) + 2j2aVaI^{0), 
showing that there is perfect separation of isospin and Ising sectors (although there is a coupling between Ising and 
spin via the term (f>(0)e). Since C3,4 involve coupling between isospin and Ising operators (/£ ^)_i and e, in this 
idealized limit one thus obtains 03^4 = identically. But in generic models, one naturally expects C3.4 ^ 1 to be finite 
but small. For this reason we ignore the operators ©3.4 in the present work. 

However, we do note that the ©3.4 operators do play a significant role in variants of the standard 2IKM which 
contain explicit coupling between isospin and Ising or spin sectors. The two-impurity Anderson model is a pertinent 
example: the low-energy effective model^iii is a 2IKM but with additional terms such as tp^{(TT^)^ ■ {Sl + Sr). These 
terms mix the spin and isospin sectors, leading thereby to indirect coupling between (/£ n)-! and e operators, and 
hence in such situations 03,4 are not expected to be small. 



B. Corrections to the Green function 



Corrections to the t matrix and conductance in the vicinity of the 2IKM critical point are directly related to 
corrections to the single particle Green functioniiS Ignoring the operators ©3^ as above, the leading correction to the 
Green function to first order in Eq. (|9]) is SG{zi, Z2) = SiG{zi, z^) -I- &2G(z\^ Z2), where 



SiG{zi,z2) = 

S2G{zi,Z2) = 



Cl 



p 



yTc Jo 

C2 



/Tr, 



dT(V'(zi)e'(0,T)V'^(z2)), 



dr(VXzi)J-i-0(O,T)V/(z2)), 



(11) 







and /? = 1/T is inverse temperature. Here we suppress the spin and channel indices, and use zi = ti + m'i, £2 = T2+ir2. 
Assuming ri > and r2 < 0, the propagator is sensitive to scattering from the impurities at the boundary located 
at r = 0. The calculation of 52G{zi, Z2) was performed by Affleck and Ludwigji^ who exploited the fact that the 
three-point function in 1^2^(2:1,22) is fully determined by conformal invariance (up to an overall constant) since J_i -cj) 
is a Virasoro primary. The electron and Ising fields, ij] and e, are chiral fields with scaling dimension 1/2, and thusiS 



52G(zi, Z2) oc 



_C2_ 



5/2 



dr- 



[sinf (zi - Z2)] 



1/2 



sin(f(T-2;i))sin(|(T-22 



3/2 ■ 



(12) 



The r integral can be expressed in terms of hypergeometric functions, which yield asymptotically the famous y/u) 
energy-dependence of the related t matrix.— 

The integrand in 5iG{zi^ Z2) is similarly a three point function, but e' ~ dr^ is not a Virasoro primary field. 
However, e itself is Virasoro primary. Pulling the derivative out of the correlation function, one obtains 



5iG{zi,z2) 



Cl 



dTa,(V(2l)e(0,T)V'^(z2)) (X 



Cl 



/ X 3/2 

[sin^(zi-Z2)]"V2(^^j drdr 



sin 



Zi) sm 



(f(^-^2) 



1/2 



(13) 
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FIG. 3: NRG energy levels SEn for odd iteration number, N. The lowest 2 levels are shown in each left/right charge and total 
spin subspace with quantum numbers {Ql,Qr, Stot) ~ (0,0,0) [dotted lines], (0,1,^) [solid lines], (1,0,^) [dashed lines] and 
(1,1,0) [dot-dashed lines]. Left panel: symmetric spin-i 2IKM with pJl ~ pJn = 0.075. Middle panel: asymmetric spin-i 
2IKM with pJl = 0.075 and pJn = 0.05. Right panel: asymmetric spin-1 2IKM with pJl ~ 0.15 and pJr = 0.05. Interimpurity 
coupling K ~ Kc tuned to its critical value in each case, and Hps = 0. 



The integrand has no singularities as long as zi and Z2 are away from the boundary, and no branch cuts associated 
with the square-root function are intersected. Periodicity in /? then implies that the integral vanishes, and so the 
contribution to the Green function from SiG{zi, Z2) also vanishes. Thus, the anomalous square-root correction to 
the Green function must come from 52G{zi, Z2) alone. In particular, our conclusion is that the coefficient of the 
square-root energy-dependence of the t matrix and hence conductance is proportional to C2, which vanishes in the 
symmetric 2IKM, obtained as Jl — > Jr. 

Higher-order corrections to the Green function can be calculated in a similar fashion. In par- 
ticular, the correction coming from the next order in perturbation theory involves integrals such as 
Jq dr dr' {'ip{zi)Oi{0,T)Oi{0,T')'ip^ {Z2)) ■ Since such integrals contains singularities when r — !• r', the correction 
does not in general vanish, even when dr and dr' are pulled outside the correlator in the case of Oi = e' . Such 
calculations are notoriously involved; and here necessitate the use of an ultraviolet cutoff 0(Tc) to avoid unphysical 
divergences. The leading correction is however expected to yield a linear energy dependence (albeit up to possible 
log corrections), in agreement with our NRG calculations, see Fig. 2 of Ref. [H 



III. 2IKM CRITICAL PHYSICS WITHIN THE NUMERICAL RENORMALIZATION GROUP 

In this appendix we consider the role of parity breaking on the RG flow, finite size spectrum and thermodynamics 
of the 2IKM. Such information can be extracted from NRG calculations)^ as now briefiy reviewed. 

Since its first application to the single-impurity single-channel Kondo model,— Wilson's NRG technique has been 
used successfully to study a wide range of quantum impurity problems (for a comprehensive recent review, see Ref. [H) . 
More recently, the increase in computing resources has permitted the detailed analysis of two-channel models involving 
several impurities, such as the 2IKM considered in the present work. 

The key element of the NRG techniqu o^^'^'^ is a logarithmic discretization of the free conduction electron Hamil- 
tonian, Hq. The continuum of states in each conduction band is divided into intervals with discretization points 
Xn = ±£)A~" (here n = 0, 1,2, ... and 2D is the bandwidth), and whose width thus decreases exponentially as the 
Fermi level is approached. A single state (the symmetric linear combination) is then retained in each interval, such 
that low-energies are exponentially-well resolved. Canonical transformation by tridiagonalization yields the Wilson 
chain representation; ^^i-'^'^ where each conduction channel corresponds to a semi-infinite chain terminated by the im- 
purity sub-system. The discretized Hamiltonian is then diagonalized iteratively: starting from the impurity, Wilson 
chain orbitals are coupled on successively and the system diagonalized. To avoid exponential growth of the Hilbert 
space, high-energy states are discarded after each step. This truncation scheme correctly allows calculation of the 
lowest-energy eigenenergies {En} of a finite Wilson chain of length N because the coupling between the Wilson chain 
orbitals N and {N + 1) scale as A~^/^. In consequence ; ^ ^1 ^ high-lying states at one iteration do not cross over and 
become low- lying states at a later iteration due to the energy scale separation inherent when A > 1. 
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FIG. 4: Uniform magnetic susceptibility TcXimpiT) vs temperature T/Tc for the 2IKM model. Shown for pJl = 0.075 > pJn, 
varying pJn = 0.075 — >■ 0.05 in steps of 0.0025, with K — Kc retuned in each case (and Hps =0), as in Fig. 2 of Ref. Hi. Thick 
solid line is the symmetric case Jl = Jr', asymmetry Jl/Jr > 1 increases in direction of the arrow. Circle points for a 2CK 
model with T}f^ = Tc- Inset shows the variation of the slope ai as a function of Jl/Jr (see Eg. I14|) . a2 ~ 0.05 is essentially 
independent of asymmetry. We have defined Tc here through T'cXimp(Tc) = 0.05. 



SE^ii) = A^/^(i?Ar(i) — min{i?Ar}) are rescaled many-particle energies (indexed i and measured with respect to the 
ground state energy), and span roughly the same energy range, independent of TV. The evolution of these levels with 
N can be understood in terms of an RG flow,— with the various fixed points giving characteristic spectra which 
do not change upon further iteration. Indeed, the lowest energy levels at the stable fixed point (obtained as TV — > oo) 
reproduce accurately the finite size spectrum obtained by CFTJ^ 

Thermodynamics can also be calculated from these NRG energy levels.— The essential step here is the identifi- 
cation of a characteristic temperature ^ A~^/^, at which thermodynamic quantities can be accurately calculated 
for a given finite iteration iV. This temperature is chosen to be high enough that the splitting of the levels incurred 
at later iterations does not affect the thermodynamic calculation; but not too high that states above the truncation 
limit contribute significantly. Thus useful physical information can be extracted from each iteration, and so the full 
temperature-dependence of thermodynamic quantities can be built up. 

A. Effect of parity-braking on RG flow and finite size spectrum 

In Fig. [3] we show the evolution of low- lying NRG eigenenergies SEn with iteration number (Wilson chain length) 
N, for 2IKMs at criticality. Specifically, we compare the channel-symmetric and -asymmetric spin-i 2IKM and the 
asymmetric spin-1 2IKM (with Hps = in each case). The RG fiow, as evidenced by the flow of these levels, is 
manifestly different for the three cases considered. In particular, levels with left/right charge and total spin quantum 
numbers (Ql, Q_R, •S't^ot) = (Oi I7 ^) ^^'^ (liO,i) are of course degenerate in the symmetric 2IKM (left panel), since 
the Hamiltonian is invariant on swapping L and R labels. By contrast no such symmetry of the bare Hamiltonian 
is present in the channel-asymmetric spin-i or spin-1 2IKMs plotted in the middle and right panels; and thus levels 
related by L o i? permutation are not in general degenerate. 

Importantly however, an emergent parity symmetry at the stable NFL fixed point is observed, with (Ql, Qr^ S^ot) 
and [Qr^Ql^S^^^) levels becoming degenerate as TV —> 00. 

Indeed, the set of NFL fixed point levels in each case is identical, demonstrating that the stable fixed point itself is 
identical, irrespective of bare model symmetries, and independent of spin-S'. The fixed point levels are indicated by 
the arrows, and correspond to the fractions 0, |, i, |, 1, as obtained for the regular symmetric 2IKM by CFTi^ 
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B. Effect of parity- braking on thermodynamics 

We turn now to thermodynamics, focusing on the 'impurity' contributio n^^i^^ to the uniform spin susceptibihty, 
XhnpC^) = (('5'tot)^)imp/r (here S^^^ refers to the spin of the entire system and (17)i,„p = (O) — {^1)q, with 
denoting a thermal average in the absence of the impurities). 

In the 2CK model, the uniform susceptibility diverges logarithmicall y at low temperatures T <C T]f^, 

n^^xl.^iT) - ai ln(r|.^^/T) + a,. (14) 

One naturally expects the uniform susceptibility of the 2IKM to behave similarly in the channel-asymmetric limit 
(with T]f^ = Tc), since here there is a mapping to the 2CK model.— However, as pointed out in Ref. [a, the uniform 
susceptibility is not singular in the regular symmetric 2IKM. Rather, it is the staggered susceptibility that is divergent 
in this case. 

In Fig.|4]we show how the uniform susceptibility TcXim-pi'r^) temperature T /Tc evolves with increasing asymmetry 
for the spin-i 2IKM. In the case of large channel asymmetry {Jl/ Jr = 1-5), the behavior is indeed that of Eq. (|14p . 
with coefficients ai and a2 essentially those of the regular 2CK model (see comparison to the pure 2CK case, circle 
points). But, in analogy to the vanishing square-root energy dependence of conductance at the symmetric point, we 
find that the coefficient ai — )• as / J/? — )• 1 (see inset). Indeed, the leading contribution to the uniform susceptibility 
in this limit can be understood^ from second-order perturbation theory in the leading irrelevant operator O2 = J_i • </>. 
From Eq. ([TU)) it then follows that ai oc {Jl — JrY for small {Jl — Jr)- Thus, there is a smooth crossover between 
the limiting cases, with divergent 2CK behavior arising in the asymmetric limit, but constant uniform susceptibility 
as T — > emerging at the symmetric point, consistent with Ref. |5|. 
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